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Abstract. We discuss several properties of static, spherically symmetric wormholes with particular 
emphasis on the behavior of causal geodesies and the propagation of linear fields. We show there 
always exist null geodesies which are trapped in a region close to the throat. Depending upon the 
detailed structure of the wormhole geometry, these trapped geodesies can be stable, unlike the 
case of the Schwarzschild black hole. We also show that test scalar fields propagating on such 
wormholes are stable. However, when a mixture of ghost and Klein-Gordon scalar fields is used as 
a source of the Einstein equations we prove that the resulting static, spherically symmetric wormhole 
configurations are linearly unstable. 
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INTRODUCTION 

Traversable wormholes has been the subject of numerous investigations over the last 
few decades. Whether they are considered as a spacetime consisting of two asymptot- 
ically flat ends connected by a throat or as a handle connecting two distinct regions of 
the same universe, wormholes exhibit properties that fascinate the public and scientists 
alike: transport of matter from one asymptotically flat end to the other and interstellar 
travel, the existence of close timelike curves and the possibility of backwards time travel, 
etc., see for example [1, 2, 3]. 

Although there exist many spacetimes describing traversable wormholes, in the con- 
text of General Relativity, a severe restriction comes from the requirement that these 
spacetimes satisfy Einstein's field equations. This leads to the necessity of exotic matter, 
ie. a configuration whose stress energy tensor violates the null energy condition [1, 4], 
and presently this type of matter is a hypothetical rather than a configuration supported 
by experimental evidences. On the other hand, the observed late time accelerated ex- 
pansion of the universe seems to demand some form of dark energy which could, in 
principle, be modeled by exotic matter, although it is not clear whether or not such 
expansion really requires matter violating the null energy condition. Therefore, so far 
neither cosmological observations nor terrestrial physics provide direct support for the 
existence of exotic matter. 

However, it has been known long ago [5] that within the framework of Quantum Field 
Theory the renormalized stress energy tensor allows the violations of all so far known 
energy conditions and this leaves open the possibility that exotic matter may be provided 
by quantum fields. These expectations have faded away after the development of the 
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quantum inequalities [6, 7, 8] that put strong constraints on the size of the region where 
violations of the energy conditions occur. For further discussions on this important issue 
and the relevance of this type of analysis for wormholes, see for example [9, 10, 11]. 
Often in the literature the existence of exotic matter is circumvented by employing 
generalized theories of relativistic gravity such as theories arising as the low energy 
limit of string theories. Examples are dilatonic Gauss-Bonnet theory or other higher 
curvature theories like f(R) theories. In these generalized theories wormholes have been 
constructed without any need for exotic matter (see for example [12, 13, 14]). 

However, irrespectively of the issue regarding the existence and nature of the exotic 
matter, an important open problem concerns the stability of wormholes. Does there exist 
a matter model in the context of General Relativity which admits stationary wormhole 
solutions which are stable with respect to small perturbations of the metric and matter 
fields? So far stability analyses have been mainly restricted to spherically symmetric 
wormholes. For the case of a massless ghost scalar field, it has been shown in [15] that 
all the static and spherically symmetric wormholes are linearly unstable, and a numerical 
evolution of the spherical nonlinear field equations reveals that these wormholes either 
expand or collapse to form a black hole [16, 17, 18]. Charged generalizations of these 
wormholes were also considered in [19] and shown to be linearly unstable. Further 
analysis includes a model involving exotic dust in combination with a radial magnetic 
field [20, 21] where linearly stable wormholes are obtained. However, as was later 
shown in [22], these configurations are unstable due to the formation of shell crossing 
singularities. Very recently, a generalization of the previous model to an exotic perfect 
fluid with pressure has been considered and claimed to yield wormholes which are stable 
with respect to radial perturbations [23]. It would be very useful to have a more general 
understanding on the stability issue of wormholes. 

The present work discusses a few properties of wormhole spacetimes which we 
believe may be of relevance to the stability problem. We start in the next section with 
an analysis of the asymptotic behavior of the causal geodesies on a given, static and 
spherically symmetric wormhole background. We show that the presence of the throat 
leads to the existence of trapped timelike and null geodesies orbiting around the throat, 
and depending upon the nature of the effective potential these geodesies can be stable in 
the sense that a small perturbation of the initial data belonging to these curves leads to 
the same trapped behavior. Next, we study the time evolution of test scalar fields on the 
wormhole background and show that they remain bounded in time. Finally, we consider 
a matter model consisting of a mixture of ghost and Klein-Gordon scalar fields, and ask 
whether or not the inclusion of the Klein-Gordon fields could act as a stabilizer agent. 
However, we find that in this model all the static, spherically symmetric wormholes are 
linearly unstable. 

GEODESIC MOTION ON A WORMHOLE BACKGROUND 

In this section we analyze geodesic motion on a given wormhole spacetime of the 
form 

M = R 2 xS 2 , g = -dt 2 + dx 2 + r(x) 2 (d$ 2 + sin 2 $ d(p 2 ) , (1) 
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where here (t,x) are Cartesian coordinates on IR 2 , (#, (p) are the standard polar coordi- 
nates on the unit 2-sphere S 2 , and r(x) is a strictly positive smooth function of x satis- 
fying limr(jc)/|jc| = 1. The function r(x) has a global minimum atx t h which describes 

the wormhole throat whose area is 4nr(x t h) 2 and which connects the two asymptotically 
flat ends at x — > +°° and x — > — °°. A particular example is r(x) = Vx 2 + b 2 , with b > 0, 
corresponding to the Bronnikov-Ellis solution [24, 25] whose properties in the context 
of wormholes were discussed in [1]. 

The causal geodesies y(A) on the spacetime (1) are described by the Lagrangian 

L = g(y,y) = -i 2 +x 2 + r(x) 2 <p 2 , 

with a dot denoting differentiation with respect to an affine parameter X. Here, we 
assume without loss of generality that the motion takes place in the equatorial plane 
■& = it 1 2. Since L does not depend explicitly on t nor on <p, the quantities E :=i and 
£ := r(x) 2 q> are constant along the trajectories. This leads to the effective mechanical 
problem 

I 2 

x 2 + V(x)=E 2 -k 2 , y(*) = _, (2) 

r[x) L 

with k 2 = for null and k 2 > for timelike geodesies. A rescaling X >->■ AX of the affine 
parameter by a positive constant A implies the transformations E h-> E/A, i >->■ £/A and 
i-)- fc/A of the conserved quantities. 



Null geodesies 



9 £ 

i+ ^2 = L ( 3 ) 



For null geodesies, the rescaling freedom mentioned above allows us to choose E = 1, 
and Eq. (2) can be rewritten as 

f 
r{xf 

Therefore, we are looking for the trajectories with energy one in the effective potential 
£ 2 /r(x) 2 which scales like the square magnitude of the angular momentum. This poten- 
tial is positive, has a global maximum at the throat, x = x t h, and decays to zero as £ 2 /x 2 
for x — > ±oo. 

In order to study the asymptotic behavior of the null geodesies, we consider an event 
p E M and a future-directed null vector 1 e T p M, and examine the asymptotic limit for 
X — > oo of the null geodesic through p with tangent 1 at p. We may parametrize 1 in terms 
of the angle a according to 

d d sin a d 

1 = + cosa— + 



dt dx r(xo) d(p ' 

where jco is the ^-coordinate of p and the normalization has been chosen such that 
i = E = 1. The condition <p = £/r(xo) 2 implies that 

£ = r(jc ) sin a. (4) 
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Before we proceed with the analysis, let us assume first that the areal function r(x) 



is strictly convex, as is the case for the Bronnikov-Ellis wormhole r{x) = \/x 2 + b 2 
mentioned above. This implies that the effective potential V(x) = £ 2 /r(x) 2 has a unique 
maximum at the throat x = x t h, see Fig. 1. With this assumption, we arrive at the 

V(x) 




Xth Xt um j n g 

FIGURE 1. The structure of the effective potential assuming that r(x) is strictly convex. The unique 
maximum corresponds to unstable circular orbits. 

following conclusions: 

(i) The radial null geodesies correspond to the choice a = or a = % for photons 
moving towards x — > +°° and x — > respectively. In this case, the angular 
momentum £ is zero, and thus the photons reach their respective asymptotic regions. 

(ii) This behavior persists for angles a close to or %, as long as the maximum of the 
effective potential, £ 2 /r(x t i 1 ) 2 , lies below one. In view of Eq. (4) this is the case if 
and only if jet j < oc c (xq) or \a — n\ < a c (xo), where 

a c (,o):=arcsin^j. 

Notice that < OC c (xq) < n/2 and <X c (xq) = %/2 if and only if the initial point is at 
the throat, xq = x t h- 

(hi) When \a\ > <x c (xq) and \a — n\ > cc c (xq), the maximum of the effective potential 
is larger than one and there is a turning point at x turn i ng such that r{x turn i ng ) 2 = £ 2 = 
r(xo) 2 sin 2 a > r(x t h) 2 . 

(iv) In the limiting case a = (X c (xq) or a = OC c (xq) + %, the maximum of the effective 
potential is exactly one and £ = ±r(x t h). If xq = x t h, the solution describes an 
unstable circular orbit, for which the corresponding light ray winds around the 
throat. When xq ^ x t h the trajectories either converge to x = ±°° or asymptotically 
approach the throat as A — > °°. 

From Eq. (3) and the definition £ = r(x) 2 <p one obtains the following expressions for the 
affine parameter X and the azimuthal angle <p: 



A(x 2 )-A(xi) 



X2 



r(x)dx 



X] 



Vr(x) 2 -£ 2 



(p(x 2 ) - <p(xi 



X2 



dx 



ix)^r(x) 2 -£ 2 



for a path starting at x\ and terminating at x%, assuming that there is no turning point 
in between. From the first expression, it follows that the affine parameter diverges as 
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%2 — > ±°° in the cases (i)-(iii) above. In the limiting case (iv), where £ 2 = r(x t / l ) 2 , 
the affine parameter also diverges when X2 approaches the throat since near the throat 
r(x) = r{x t h) + jr" (x t h)(x — x t h) 2 + G(x — x t h) 3 . Therefore, the effective potential is 
complete [26] in the sense that each orbit exists for infinite values of the affine parameter. 
In particular, this implies that the wormhole spacetime (1) is null geodesic complete. 

We may summarize the results so far in the following way: Let us introduce the 
effective phase space r := E x S 1 of possible initial data (jco, e M ) for the null geodesies. 
Consider the subset 

r„:={(xo t (**)er: lim \x(X)\ = °°} 

of those data which give rise to a null geodesies which escapes to infinity in the future. 
According to the analysis above, the complement of this set, describing future trapped 
null geodesies, is characterized by 

r trapped :=r\r^ = {(xo^)er:a=l ( a ^l 5°™^** U 

pp v \ v ' \a c (x )+n for xq> xth jj 

Since this set has measure zero in T, it follows that Trapped is unstable in the sense that 
a generic perturbation of initial data (xo, a) in this set leads to an orbit which reaches 
infinity in the future. 

The results so far are based on the assumption that the areal function r(x) is strictly 
convex, implying that the effective potential has a unique local maximum at the throat. 
However, consider now a background wormhole for which the effective potential ex- 
hibits the structure shown in Fig. 2. Such potentials arise in wormhole models [27] that 
were proposed as black hole foils [28]. In contrast to the previous case where the ef- 



V(x) 




Xi x 2 



FIGURE 2. An example for an effective potential with two local maxima and one local minimum, 
corresponding, respectively, to unstable and stable circular orbits. There is a potential well between X\ 
and X2 giving rise to oscillatory motion. 

fective potential has a unique local maximum, in the present case there is a potential 
well between x\ and X2 with a local minimum corresponding to stable circular orbits. 
Let us consider initial data (xo,e' a ) lying inside the well, that is x\ < xq < xj. It follows 
using similar arguments than in the previous case that for initial angles a satisfying 
\a\ > (Xi(xq) and \a — iz\ > (Xi(xo) (ie. those angles sufficiently close to n/2 or 3n/2) 
with 



CCi (xq) = arc sin 
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the motion is oscillatory, and the corresponding light rays follow a helicoidal type 
trajectory. Since these trajectories are trapped, the set 

| (x , e ia ) e r : x\ < xq < x 2 , | sin a | > ^—^ j c T trapped 

is an open subset of T trapped- As a consequence, = T \ T trappec i cannot be dense in T, 
and it is not true that a generic null geodesic escapes to infinity, like in the previous case. 



Timelike geodesies 

For timelike geodesies it is natural to choose the scaling freedom described below 
Eq. (2) such that k = 1 , implying that the affine parameter measures proper time. The 
initial four-velocity u can be parametrized in terms of an hyperbolic angle % > and a 
conventional angle a according to 

d . , / d sin a d 
u = cosh^— + sinh^ cos a— + 



dt \ dx r(xo) d(p 

where xq is the x-coordinate of the initial point p. The conservation laws i = E and 
r(x) 2 (p = i imply that 

£' = cosh^, £ = r(xo) smashing, (5) 
and consequently, the geodesic equation (2) can be rewritten as 

x 2 +(^7^) sinh 2 ^sin 2 a = sinh 2 ^. (6) 
V r{x) J 

Therefore, concerning the asymptotic limit of the timelike geodesies, we may draw 
similar conclusions than in the case of null geodesies. Now the effective phase space 
describing the set of possible initial data (xo 7 e x 7 e ia ) is T = R x (1,°°) x S , and we may 
write it as the disjoint union T = T trappec i OToo, were 

T ao :={(xo,e x ,e ia )eT: lim |x(A)|=oo} 

and Ttrapped '■= T \ Too, as before. For an effective potential of the form displayed in 
Fig. 1 it follows again that T trappe d is unstable, implying that there are no stable trapped 
timelike geodesies. In contrast to this, for the example shown in Fig. 2 the set T tra pped 
contains the open subset 

S^(xo,e x ,e ia ) eT-.xi <x <x 2 ,X >0,|sina| > ^yj, 
describing stable trapped timelike geodesies (planetary motion). 
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EVOLUTION OF FIELDS ON A WORMHOLE BACKGROUND 



After having analyzed geodesic motion on the wormhole background (1), describing 
the evolution of photons or massive test particles, in this section we make a few remarks 
about the propagation of test fields on the background spacetime (1). For definiteness, 
we consider here a scalar field 4> of mass /i whose dynamics is described by the Klein- 
Gordon equation 



-V%4> + jU 2 4> = 0. 



(7) 



When specialized to the spherically symmetric background (1), this can be reduced to 
a family of radial equations through a decomposition 4> = r(x)~ l £ (j)e m (t,x)Y em (&, <p) 



into spherical harmonics Y 



im 



d 2 <l>£m 
dt 2 



r"{x) £(£+\) 2 
r[x) r[x) L 



(8) 



The first term, r"(x)/r(x), in the effective potential Vm comes from the curvature of 
the spacetime background (notice that this term would be zero in the Minkowski case 
r(x) = x), the second term is the centrifugal term and is the quantum analog of the 
effective potential in Eq. (2) for the geodesic motion, while the last term is just inherited 
from the Klein-Gordon equation (7). 

For the Bronnikov-Ellis wormhole r(x) = \Jx 2 + b 2 , for instance, the curvature term is 
r"{x) I r(x) = b 2 / r(x) 4 and Vg^ is positive and has a global maximum at the throat x = 0. 
Therefore, it describes a potential barrier and in principle, one can estimate transmission 
and reflection coefficients for scattering processes involving incoming radiation from 
one asymptotic flat end. It would be interesting to perform such an analysis and to 
compute the quasinormal modes which are expected to arise from the poles of the 
analytic continuation of these coefficients [29] . 

Solutions of Eq. (7) generated by smooth initial data with compact support remain 
bounded in time. This can be inferred from the existence of the conserved energy 



RxS 2 



-^-J +f j Di<Z>Dj<f> + iJ, 2 <S> 2 



Vfd 3 x, 



(9) 



where y = dx 2 + r(x) 2 (di} 2 + sin 2 # d<p 2 ) is the three-metric induced from g on the 
spatial, hypersurface-orthogonal slices t = const. This energy is the conserved quantity 
on such time slices, belonging to the conserved four-current 7^ = —T^ v k v which is 
constructed from the stress-energy tensor T^ v of the scalar field and the static Killing 
vector field k = d /dt of the spacetime (1). Due to the orthonormality of the spherical 
harmonics, this energy reduces to a sum E = £ E( m with 



oo - 




2 

+ r{xf 




f §tm\ 


/ 

-co L 


dt 


-( 
dx ' 


\r(x)J 



+ 



1(1+1) 

r{x) 2 



<?m| 2 +M 2 |<fen| 2 



dx, (10) 
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when decomposed into spherical harmonics. Since each term Ei m is positive definite, 
one obtains a bound on the L 2 norm of §{ m and their first order derivatives, and from this 
a uniform bound on <$>£ m } 

In fact, similar conclusions can be obtained for Maxwell's equations on the wormhole 
background (1). It follows that the time evolution of scalar and electromagnetic fields 
from smooth initial data with compact support on such backgrounds stays bounded, 
and in this sense these fields are stable. However, as discussed in the next section, the 
situation may change drastically when the self-gravity of the fields is taken into account. 

WORMHOLES SUPPORTED BY SCALAR FIELDS AND THEIR 

STABILITY 

We consider here a family $ = A = 1,2, ... ,ra, of self-gravitating, noninteract- 
ing, minimally coupled, massless scalar fields. The action is 

S[%M = J (-T^ + ^ V ^ A ' V M^) V=gd\ (11) 

where R and V denote the Ricci scalar and the covariant derivate, respectively, 
associated to the spacetime metric g, G is Newton's constant and (&ab) = 
diag(— 1,— 1,..., — 1, + 1,..., + 1) is a flat metric on internal space with r negative 
and s positive entries, corresponding to r phantom and s Klein-Gordon scalar fields. The 
corresponding field equations are 

R^v = %KGh AB V^ A -V v f, (12) 
= V%0 A , (13) 

where R^ v is the Ricci tensor associated to g. The action (11) and field equations (12,13) 
are invariant with respect to the global transformations 

(j) A ^A A B ({) B + a A , AeO(r,s), a G R m , (14) 

in the internal space. Using this symmetry, particular solutions of Eqs. (12,13) can be 
constructed by applying the transformation (14) to a solution having only one component 
in internal space, 4> = (0, 0, 0, . . . , 0), in which case the system reduces to the Einstein- 
ghost scalar field equations whose solutions in the static, spherically symmetric sector 
are well known [24, 25]. In fact, we claim that all the static, spherically symmetric 
wormhole solutions of Eqs. (12,13) can be obtained by this method. In order to show 
this, assume (g, <J>) is such a solution. We can write the metric in the form 

g = - e -M*) dt 2 + e 2a( X ) dx 2 + ^2 ^2 + ^2 # ^2) ? 



This follows from a classical Sobolev estimate, see for instance [30]. With a little bit more work one 
could probably also show that 4> itself is uniformly bounded on the spacetime (1) provided the initial data 
is smooth and has compact support. 
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with a regular function a(x) and r(x) the areal function. In this parametrization, the wave 
equations (13) reduce to 



d_ 
dx 



r(x) 2 e' 



dx 

whose solutions are described by 



0, A=l,2,...ra, 



f e 2a ^ 

A (x)=/(x)a A + j8 A , f(x)= J -^dy, A = l,2,...m, 

— oo 

with integration constants a A and f5 A . The important point to notice here, is that the 
function f(x) is common to all the A 's. Therefore, we can apply the symmetry (14) 
to transform (j) A to a vector in internal space of the form 4>' — ((j)', 0, 0, . . . , 0), provided 
we ensure that the vector (a A ) is timelike in internal space, hAB<X A <X B < 0. In order to 
show that this must be the case, we consider the combination R — 2R X X of the Einstein 
equations (12) which yields 

r\x) 2 , V(jc) e 2a W (d$ A \ (d<$> 



r{xf 



Evaluating at the throat x = x t h and using r '(x t h) = combined with the regularity 
of a(x), it follows that AnGh AB a A a B f (x th ) 2 = -e 2a( - x ' h ^ /r(x th ) 2 , which shows that 
f'(xth) 7^ an d that hp, B a A a B is negative. Therefore, as a consequence of the symme- 
try (14), any static, spherically symmetric wormhole solution of Eqs. (12,13) is equiva- 
lent to a static, spherically symmetric solution of the Einstein-ghost scalar field theory. 
Since the latter are known to be linearly unstable [15] it follows that any static, spheri- 
cally symmetric wormhole solution of the theory involving r ghost and s Klein-Gorden 
scalar fields is unstable with respect to linear fluctuations of the metric and scalar fields. 
The numerical results for the nonlinear evolution in [16, 17, 18] show that such worm- 
holes may collapse to a Schwarzschild black hole. 

In conclusion, this model shows that the inclusion of one or more Klein-Gordon fields 
to the system does not lead to wormholes which are linearly stable. This result, in 
combination with the effect of the electric charge on the stability [19], suggests that 
it might be unlikely to stabilize a wormhole by adding ordinary matter. 



CONCLUSIONS 

In this work we have presented some basic properties of static, spherically symmetric 
wormhole spacetimes. For the wormhole models given in Eq. (1) we have shown null 
and timelike geodesic completeness. Furthermore, we have shown that depending upon 
the structure of the effective potential, there might exist trapped timelike and null 
geodesies winding around the throat. In contrast to the Schwarzschild black hole, cases 
are exhibited involving stable trapped null geodesies. As far as fields are concerned, we 
have shown that a test scalar field on the wormhole background (1) remains bounded 
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in time for smooth initial data with compact support. Finally, we consider a mixture of 
self-gravitating ghost and Klein-Gordon scalar fields and prove that the resulting static, 
spherically symmetric wormholes are linearly unstable. In this context, it is interesting to 
mention that the nonlinear stability of Minkowski spacetime in the Einstein-scalar field 
model has been established not only for a Klein-Gordon field but also for a ghost scalar 
field (see the comments and references in appendix B.5 of [31]). Therefore, it seems 
that both the self-gravity and the topology are key factors leading to the instability of the 
wormholes. 
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